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Light gauge bosons can lead to resonant interactions between high-energy astrophysical neutrinos 
and the cosmic neutrino background. We stndy this possibility in detail, considering the ability 
of IceCube to probe such scenarios. We find the most dramatic effects in models with a very 
light Z' {mz' ^ 10 MeV), which can induce a significant absorption featnre at ~ 5-10 TeV x 
(mz'/MeV)'^. In the case of the inverted hierarchy and a small snm of neutrino masses, snch a 
light Z' can result in a broad and deep spectral featnre at ~0.1-10PeV x (m^z/MeV)^. Current 
IceCube data already excludes this case for a Z' lighter than a few MeV and conplings greater than 
g ~ 10“"^. We emphasize that the ratio of nentrino flavors observed by IceCube can be nsed to 
further increase their sensitivity to Z' models and to other exotic physics scenarios. 


I. INTRODUCTION 


In addition to providing a new window into the ori¬ 
gin of the cosmic ray spectrum, the observation of as¬ 
trophysical neutrinos allows us to probe fundamental 
physics. More specifically, IceCube’s recent detection of 
high-energy astrophysical neutrinos enables us to study 
and constrain a range of phenomena at higher energies 
and over longer baselines than can currently be tested 
in laboratory environments. The flavor ratios of high- 
energy astrophysical neutrinos can be used to constrain a 
variety of new phenomena, including neutrino decay and 
Lorentz violation [1-18]. Alternatively, a few PeV neu¬ 
trino interacting with a nucleon at rest has a center-of- 
mass energy of a few TeV, enabling high energy neutrino 
telescopes to constrain a range of TeV-scale physics sce¬ 
narios, including TeV-scale gravity models [19-23], and 
models featuring TeV-scale leptoquarks [24, 25]. These 
observations also make it possible to study interactions 
between high-energy neutrinos and the cosmic ray neu¬ 
trino background (Ci^B). In particular, scenarios featur¬ 
ing a light gauge boson have received some recent atten¬ 
tion within this context [26-34]. 

A new gauge boson with couplings to Standard Model 
neutrinos will induce a scattering resonance with the CvH 
at an energy given by: 
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For even very small couplings, such a resonance can lead 
to the efficient absorption of high-energy neutrinos over 
cosmological distances. 

In this paper, we revisit the possibility of using Ice- 
Cube (or future high-energy neutrino telescopes) to 
search for the effects of an MeV-scale gauge boson on 
the high-energy cosmic neutrino spectrum. In doing so, 
we consider the impact on both the shape of the neutrino 
spectrum, as well as on the ratio of flavors that reach the 
Earth. In the window of parameter space that is capable 


of explaining the measured value of the muon’s anoma¬ 
lous magnetic moment, significant effects can result from 
such a Z'. 

In the following two sections, we review IceCube’s dis¬ 
covery of high-energy astrophysical neutrinos, and sum¬ 
marize the motivations for a model with an MeV-scale Z' 
with couplings to Standard Model neutrinos. In Sec. IV, 
we describe the interactions mediated by such a Z' be¬ 
tween high-energy neutrinos and the cosmic neutrino 
background. The impact of such interactions on the spec¬ 
trum and the flavor ratios of the high-energy astrophys¬ 
ical neutrino flux is discussed in Sec. V. In Sec. VI, we 
summarize our results and conclusions. 


II. ICECUBE’S OBSERVATION OF 
HIGH-ENERGY ASTROPHYSICAL NEUTRINOS 

Recently, the IceCube Collaboration has reported the 
observation of a diffuse flux of high-energy extraterres¬ 
trial neutrinos, consisting of 37 neutrino candidate events 
with energies ranging from 30 TeV to 2 PeV [35-37]. 
Although the origin of these neutrinos is currently un¬ 
known, they appear to be approximately isotropically 
distributed across the sky [37], suggesting an extragalac- 
tic origin (see, however, Refs. [38-44]). The spectrum of 
these particles is well fit by a power-law with an index of 
7 = -2.6 [45]. 

Several features of this neutrino population are sug¬ 
gestive of a connection with the cosmic ray spectrum. 
In particular, the generation of the observed neutrino 
flux requires that ~20% of the PeV-EeV protons accel¬ 
erated by cosmic ray sources undergo photo-meson in¬ 
teractions; a fraction that could easily be accommodated 
in realistic astrophysical environments [46]. Stated an¬ 
other way, the observed neutrino flux is below, but not 
very far below, what is known as the “Waxman-Bahcall 
bound” [47, 48]. Furthermore, the numbers of showers 
and muon track events observed at IceCube is consistent 
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with a flavor ratio oi Uf, '■ Vfi ■ Vt = ^ ^ ^ (although with 

large error bars) [45, 49], consistent with that predicted 
from photo-meson interactions after accounting for os¬ 
cillations. Several plausible classes of sources have been 
proposed for these neutrinos, including active galactic 
nuclei [46, 50-53], starburst or star-forming galaxies [54- 
57], and low-luminosity gamma-ray bursts [46, 58-60]. In 
addition to these more conventional astrophysical source 
classes, exotic origins for IceCube’s neutrinos have also 
been considered, such as the decays or annihilations of 
long-lived superheavy particles [61-70]. In this study, we 
remain agnostic as to the specific origin of these neutri¬ 
nos, assuming only that they originate from extragalactic 
sources. 


III. A LIGHT Z' WITH COUPLINGS TO 
NEUTRINOS 


New gauge bosons appear within many new physics 
scenarios [71]. For example, additional broken Abelian 
1/(1) gauge symmetries and the Z' bosons that accom¬ 
pany them are predicted by many Grand Unified Theo¬ 
ries (GUTs), including those based on the groups S'O(IO) 
and Eq [72, 73]. New massive gauge bosons also appear 
within the context of many string inspired models [74- 
82], little Higgs theories [83-86], dynamical symmetry 
breaking scenarios [87-89], models with extra spatial di¬ 
mensions [90-93], and many other proposed extensions 
of the Standard Model [94-96]. 

In this paper, we are primarily interested in light gauge 
bosons {mz> < 1 GeV) with nonzero couplings to Stan¬ 
dard Model neutrinos. Gonstraints on such a particle’s 
couplings to electron neutrinos are quite stringent, how¬ 
ever, motivating us to focus on models in which the Z' 
couples only to 2nd and/or 3rd generation leptons. 

If the couplings of a Z' are assigned arbitrarily, anoma¬ 
lies are generally introduced, violating the principle of 
gauge invariance. To construct a self-consistent theory, 
care must be taken to ensure that all such anomalies 
cancel. A well-known example of an anomaly-free Z' is 
that arising from the gauge group This is the 

only anomaly-free 1/(1) group with nonzero charge as¬ 
signments to Standard Model neutrinos that can lead to 
an experimentally viable MeV-scale Z' without requir¬ 
ing the addition of any exotic fermions. A 1/(1) group 
charged under only muon or tau number is also a possi¬ 
bility, although new chiral fermions must be introduced 
in these cases, charged under both SU{2)\y and t/(l)Y, 
as well as under the new U{1)^ or U{\)r group. 

A light Z' with couplings to the muon can be mo¬ 
tivated by the measurement of the muon’s anomalous 
magnetic moment, which currently differs from the value 
predicted by the Standard Model with a significance of 
approximately 3.6 ct [97]. With efforts currently underway 
to improve this measurement [98, 99] and to reduce the 
related theoretical uncertainties [100-105], it should be¬ 
come clear within the next several years whether or not 


this is an authentic sign of new physics. Among other 
possibilities (see, for example. Refs. [106-111]), an MeV- 
scale Z' with small couplings to the muon could plausibly 
account for this measurement. 

A Z' with a vector coupling to muons, g^, leads to 
the following contribution to the muon’s magnetic mo¬ 
ment [112, 113]: 
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For mz' ^ the measured value can be accommo¬ 
dated for ~ (3 — 6) X 10“^, whereas for mz' —1 GeV, 
couplings an order of magnitude larger are required. Al¬ 
though this parameter space is in conflict with measure¬ 
ments of muon pair production in muon neutrino-nucleus 
scattering for mz' ^ 500 MeV [27, 114], lower values 
of mz' remain viable (below 1 MeV, constraints from 
big bang nuclosynthesis and the cosmic microwave back¬ 
ground can also be relevant [28, 115-117]). 

As stated above, a Z' resulting from the U{1)^ or 
U(l)r groups requires the introduction of new chiral 
fermions charged under SU{2)w and U(\)y in order to 
cancel anomalies. Furthermore, the masses of these ex¬ 
otics are bounded by the requirement of perturbativity, 
which requires [118]: 


^exotic ^ 108 GeV X 


mz' \ 
lOMeVy 


0.0005 \ 
9z' ) 



(3) 


where Zip is the charge assignment for the scalar whose 
VEV breaks the U{1) responsible for the Z'. From this, 
we learn that the required exotics must be rather light, 
and will be subject to constraints from accelerators. 


IV. HIGH ENERGY NEUTRINO 
INTERACTIONS WITH THE COSMIC 
NEUTRINO BACKGROUND 

In generality, a Z' boson will couple to Standard Model 
neutrinos through the following interaction: 

Cz'v = gz'QagZ'^Va-f^^PLVg (4) 

= gz'QijZ'^Vijf^PLi^j, (5) 

where neutrinos with Greek (Latin) indices refer to the 
flavor (mass) basis. The C/(l)' gauge coupling is gz' and 
the charge assignments are contained in the matrix, Qap- 
For example, the charge matrix for [/(l)^_r is given by 
Qap = diag(0,1,-1). In the mass basis, the charge ma¬ 
trix is represented as = Ul^^QapUgj. Where Uai is 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix 
which rotates from the mass to the flavor basis. Here we 
have assumed that the charged leptons do not mix, such 
that the neutrino mixing matrix is entirely determined by 
the PMNS matrix. While this is a convention in the Stan¬ 
dard Model, deviating from this assumption does lead to 
observable effects in models with gauged lepton numbers 
[119]. 
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This interaction results in Vi — Vj scattering, with the 
following cross section: 


The thermal averaging in the second term of Eqn. 7, 
and analogously in the third term, is given by: 


a{viVj 


kl 


3tt (s — m^/)^ + w^/T^; ’ 


( 6 ) 

where the sum is over the final neutrino states, k and £. 
For example, = 2 for the 17(1)^-,- case. The t- 

channel contribution to the cross section is ignored, as it 
is highly suppressed relative to the s-channel resonance. 
The width of the Z' into neutrinos is given hy Fz' = 
J2ij Qij5l'™z'/247r. In the parameter space of interest 
to IceCube, the mass of the Z' is less than 2m^, allowing 
us to safely neglect decays to muons or taus. 

To calculate the spectrum of neutrinos at Earth, we 
solve the following coupled set of integro-differential 
equations [28, 29, 120]: 


,, .Hiz) drii ^ ^ , 

— (1 + z) - = JiiEg, z) 

c dz 
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Here, El{z) is the Hubble parameter as a function of red- 
shift, Eq is the neutrino energy as measured at Earth, 
and ni,{z) is the proper number density of neutrinos in 
the CvB. The hrst term on the right-hand side accounts 
for the source spectral and density evolution through cos¬ 
mology. The second term accounts for neutrinos of state 
i, scattering off a thermal distribution of CvB neutrinos 
of state j, thereby attenuating the neutrino flux. The 
last term accounts for the regeneration of the scattering 
products from the process described in the second term. 
Here a neutrino of state k, with energy E', scatters with 
the CvB to produce two neutrinos, one of which is of 
state i with energy Eq. Since this process is s-channel, 
the differential cross section can be broken into the total 
cross section and a distribution function in the outgoing 
neutrino energy space: 

^iE',z)=akj{E',z)f{E',Eo), (8) 

arjQ 

where 
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Note that this differential cross section accounts for both 
outgoing neutrinos, such that / dEof{E',Eo) = 2. Ex¬ 
pressing the differential cross section in this way simpli¬ 
fies the numerical implementation of Eqn. 7. 
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where the momentum dependence of the scattering cross 
section can be accounted for by evaluating Eqn. 6 with 
s = 2 Eq{ 1 + z){\/m'j -\-p^ — \p\ cos6l). In the limit that 
the neutrinos are much heavier than the effective temper¬ 
ature, the usual form is recovered. We take the effective 
temperature of the CvB at z = 0 to be Tq = 1.7 x 10“^ 
eV. 

While an absorption resonance is predicted for each 
neutrino mass eigenstate, very light neutrinos (with 
masses comparable to or less than the effective temper¬ 
ature the thermal distribution of the CvB) will lead to 
broad spectral features. In scenarios in which the lightest 
neutrino is nearly massless, the energy of the resonance 
is determined by the energy (and therefore the tempera¬ 
ture) of the neutrino, rather than by its mass. One then 
expects scattering to be most relevant at ~ m\,/2TQ} 

Throughout this paper, we adopt values for the neu¬ 
trino mass splittings and mixing angles as presented in 
Ref. [121], with the exception of the CP violating an¬ 
gle which we take to be i5=0 (which is well within 2 ct of 
the central value). For the cosmological parameters iJo, 
Firm s-nd Fl\, and for the upper limit on the sum of the 
neutrino masses, we adopt the values presented by the 
Planck Collaboration [122]. 


V. THE IMPACT OF A LIGHT Z' ON THE 
SPECTRUM AND FLAVOR RATIOS OF THE 
HIGH-ENERGY NEUTRINO FLUX 

A. Main Results 

In this section, we present the results of our calcula¬ 
tions. As our primary model of interest, we consider a 
Z' associated with the gauge group [/(l)p_r. For each 
value of mz', the coupling gz' is chosen to match the 
measured value of the muon’s magnetic moment, Aa^: 


4 X 10-4 

for mz> 

= 1 MeV 


5 X 10-4 

for mz> 

= 10 MeV 

(10) 

8 X 10-4 

for mz> 

= 100 MeV. 



For simplicity, we hrst consider the case in which all 
of the neutrino sources are located at a redshift of z = 1, 


^ More precisely, one expects the cross section to peak at Ei, ^ 
m|,/2(p>, with (p) = Si 3.15r. 
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Normal, Smv=i0.058 eV 
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mz =10 MeV 
mz.=100 MeV 
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Inverted, Lmv=i0.10 eV 
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FIG. 1. The spectrum of neutrinos at Earth, after including the effects of a Z' associated with the gauge group 
The coupling of the Z' has been chosen in each case to accommodate the measured value of the muon’s anomalous magnetic 
moment {gz' = 4 x 10“"*, 5 x 10“^^ and 8 x 10“^ for m^/ = 1, 10 and 100 MeV, respectively). Here, we have assumed a 
population of sources at 2 = 1 which inject neutrinos with a power-law spectrum of index of -2.6, and with an initial flavor 
ratio of i/e : = 1 : 2 : 0 (as predicted from pion decay). We show results for the normal and inverted hierarchies, and for 

maximal and minimal values of the sum of the neutrino masses. 


and which emit a spectra characterized by a power-law 
of index 7 = — 2 . 6 : 

J{E,,z) = AE2Siz-l), (11) 

where the normalization, A, is chosen to fit the IceCube 
data. We adopt an initial flavor ratio of ■. ■. 1 /^ = 

1 : 2:0 (as predicted from pion decay), which rapidly 
evolves to approximately 1 : 1:1 via oscillations. 

In Fig. 1, we plot the neutrino spectrum at Earth, in¬ 
cluding the effects oi a. Z'. Results are shown for three 
choices of mz', and for both normal and inverted hi¬ 
erarchies, as well as maximal and minimal values for 
the sum of the neutrino masses. For each hierarchy, we 
take 'Y2171^=0.23 eV to be the maximal value allowed by 
cosmological constraints [122]. For the minimal sum of 
masses, we adopt 0.058 eV and 0.10 eV for the normal 
and inverted hierarchies, respectively. 

In each case, the neutrino spectrum is altered by the 
interactions of the Z', although in ways that vary consid¬ 
erably depending on the mass of the Z' and on the hier¬ 


archy and the masses of the neutrinos. In the case of the 
normal hierarchy with ^ m,y=0.23 eV, for example, an 
absorption feature appears at E,, ~ 5 TeV x {mz' /MeV)^. 
This feature results from the scattering with all three 
neutrino mass eigenstates, and the individual resonances 
cannot be easily distinguished. In the normal hierar¬ 
chy with the minimal sum of masses, two features are 
visible: one at Ei, ~ 10 TeV x (ruz'/MeV)^ from scat¬ 
tering with the heaviest mass eigenstate, and another at 
El, ~ 50 TeV X (m^'/MeV)^ resulting from the combi¬ 
nation of the two lighter mass eigenstates. 

The results are somewhat different in the case of 
the inverted hierarchy. As found for the normal hi¬ 
erarchy, an absorption feature appears at Ei, ~ 5 — 
10 TeV X (m^'/MeV)^. In the case of ^to,,= 0.23 eV, 
this is the collective consequence of all three mass eigen¬ 
states. For ^TOi/=0.10 eV, however, this feature is 
induced only through scattering with the heaviest two 
eigenstates. In this case, the lightest neutrino leads in¬ 
stead to a very broad and potentially deep spectral fea- 
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FIG. 2. The fractions of neutrinos of each flavor at Earth, including the effects of a Z'. The models and other assumptions are 
the same as adopted in Fig. 1. 


ture, covering a range of energies between ~ 0.1- 
10 PeV X {mz' /MeV)'^. This broad absorption feature 
is most clearly visible in the case of mz' = 1 MeV. In 
this case, the spectrum reported by the IceCube Collab¬ 
oration (shown as error bars) is incompatible with a Z' 
lighter than a few MeV. Note that this feature does not 
appear in the case of the normal hierarchy because the 
lightest neutrino is largely of i/e flavor, and thus does 
not couple to the Z' under consideration. In contrast, 
the lightest mass eigenstate in the inverted hierarchy is 
primarily composed of and Vt- 

In addition to the impact on the high-energy neutrino 
spectrum, a light Z' can alter the ratio of neutrino fla¬ 


vors that reach Earth. In Fig. 2, we plot these ratios for 
the same range of scenarios considered in Fig. 1. Simi¬ 
lar to the spectrum, the most significant effects are seen 
in the case of the inverted hierarchy with a sum of neu¬ 
trino masses near the minimal value. In this case, the 
relative flux of electron (muon) neutrinos is suppressed 
(enhanced) over a wide range of energies, especially for 
the case of < 10 MeV. Such an extreme departure 
from astrophysical expectations could plausibly be tested 
in the future by IceCube. 

Constraints on the flavor ratios of IceCube’s neutrinos 
have been placed by comparing the distribution of muon 
track events (generated in charged current interactions 
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FIG. 3. A comparison of the neutrino spectra predicted for 
a Z' associated with the U{l)fi,-T (blue), 17(1)^ (red), and 
U{1)t gauge groups. Results are shown for the case of an 
inverted hierarchy, eV, and mz' = l MeV. 

of and i?^) to the distribution of showers (generated 
most efficiently in charged current interactions of Ve, i’ej 
Vt and Dr, as well as in neutral current interactions of 
all flavors). The results of this comparison have thus far 
been compatible with a ratio oi ^ ^ 

although with large error bars [45]. More specifically, 
extreme ratios of r'e : : i/r = 0 : 1 : 0 and Ve '■ v^j, '■ Vr = 

1:0:0 have been excluded at the level of S.Str and 2.3cr 
significance, respectively. 

Future flavor ratio measurements by IceCube are ex¬ 
pected to be strengthened by improvements in veto tech¬ 
niques, and by searches for events unique to tau neutri¬ 
nos.^ Showers generated via the Glashow resonance (at 
= m^/2me « 6.3 TeV) also provide an opportunity 
to constrain the electron antineutrino fraction at very 
high energies. 



FIG. 4. The fractions of neutrinos of each flavor predicted 
for a Z' associated with the t/(l)^ (upper frame) and U{1)t 
(lower frame) gauge groups. Results are shown for the case 
of an inverted hierarchy, eV, and m^'^l MeV. 


B. Alternative Gauge Groups 

Thus far, we have restricted our calculations to the 
case of a Z' associated with the gauge group 17(l)^_r. 
As discussed in Sec. Ill, however, one could also con¬ 
sider mediators arising from other gauge groups, such 
as U{l)fi or U{l)r- In Fig. 3 we compare the neutrino 
spectrum predicted for each of these three choices, for 
the case of the inverted hierarchy with ^m,y«0.I0 eV 
and 771^'=I MeV. From this comparison, we find that 


the U{l)fj,-r model results in a somewhat smaller degree 
of absorption at very high energies, despite the fact that 
more neutrino flavors participate in the scattering. This 
somewhat counterintuitive result is due to a cancellation 
between the and r contributions to < 533 . In Fig 4, we 
plot the flavor ratios in each of these scenarios. Here, the 
magnitude of the impact of the Z' is lessened relative to 
that predicted in the U{l)^-r case (see Fig. 2). 


C. Source Distributious 


^ At energies above a few PeV, tau neutrinos and antineutrinos can 
generate a tau lepton that travels an observable distance before 
decaying; the mean distance traveled is Rr = ErCTr jrrir ~ 50 m 
xEr/flPeV). As a result, very high-energy tau neutrinos can 
yield events with two showers (double bang events) [123, 124], as 
well as events with one observed shower followed by or preceded 
by a tau-induced track (lollipop events) [2]. 


Up to this point, our calculations have taken all of 
the neutrino sources to reside at a distance of z = 1. 
This distribution was adopted for simplicity, and reflects 
a plausible average for the distance traveled by a neu¬ 
trino detected by IceCube. In this subsection, we con¬ 
sider more realistic redshift distributions for the sources 
of the high-energy neutrinos observed by IceCube. 

The first possibility we consider is a distribution with 
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FIG. 5. A comparison of the neutrino spectra predicted for 
four source distribution models: all sources at 2 = 1 (dashed 
blue), sources with a constant comoving number density out 
to 2inax = 1 (short dashed red), 2max = 3 (dot-dashed purple), 
and a distribution of sources that follows the star formation 
rate (solid orange). Results are shown for the case of an in¬ 
verted hierarchy, 5])m^Ri0.10 eV, and mz/ = l MeV. 

a constant comoving number density out to z^ax) beyond 
which no sources exist: 

JcomovingiE^, z) = A {1 + zf 0(Zmax-2:), (12) 

where again we take 7 = —2.6. The second distribution 
that we consider follows the star formation rate [125]: 




Jsfr{E.,z) = AE2il + zf fsFRiz), (13) 

where 

({1 + z)^-^ z<l 

fsFRiz) = I 23-’ (1 + z)-0-3 1 < z < 4 (14) 

1^23.7 53.2 (i^^)-3.5 ^> 4 _ 

For sources distributed according to the star formation 
rate, we consider redshifts up to z = 7, beyond which 
contributions to the final spectra are negligible. The neu¬ 
trino spectra predicted from these distributions are com¬ 
pared in Fig. 5, as are their flavor ratios in Fig. 6 . From 
these figures, we see that our results are qualitatively 
insensitive to the precise choice of redshift evolution. 
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FIG. 6. The fraction of neutrinos of each flavor predicted 
for three source distribution models: sources with a constant 
comoving number density out to Zmax = 1 (upper frame), 
Zmax = 3 (middle frame), and a distribution of sources that 
follows the star formation rate (lower frame). Results are 
shown for the case of an inverted hierarchy, eV, 

and mz/=l MeV. 




















VI. SUMMARY AND CONCLUSIONS 

IceCube’s recent detection of high-energy astrophysi- 
cal neutrinos provides us with an opportunity to study 
the interactions of these particles at higher energies and 
over longer baselines than are currently possible in labo¬ 
ratory environments. In this paper, we have considered 
how light (~I-100 MeV) gauge bosons, with couplings 
to Standard Model neutrinos, could impact the spectrum 
and flavors of the neutrinos observed by IceCube. 

New gauge bosons are predicted within a variety of ex¬ 
tensions of the Standard Model. Of particular interest is 
the Z' that arises from the anomaly-free 17(1 )^_t- gauge 
group. For masses in the range of mz' ~ 1 — 500 MeV 
and a coupling of gz' ~ 10“^, such a particle can explain 
the measured value of the muon’s anomalous magnetic 
moment, without conflicting with constraints from accel¬ 
erators or cosmology. For this range of masses and cou¬ 
plings, high-energy astrophysical neutrinos can scatter 
resonantly with the cosmic neutrino background, lead¬ 
ing to absorption features in the spectrum observed at 
Earth. By measuring the spectrum and the flavor ratios 
of the extragalactic neutrino flux, IceCube can constrain 
or provide evidence for such models. 

We found the most dramatic effects in models with 
a very light Z' {niz' ^ 10 MeV), which induces a 
significant absorption feature at Ej, ~ 5 — 10 TeV x 
(m^'/MeV)^. Although not currently constrained, such 
a feature could plausibly be measured by IceCube or by 
next generation neutrino telescopes. Furthermore, in the 
case of the inverted hierarchy with the lightest neutrino 
lighter than ~10“^ eV, such a Z' can lead to a very broad 


and deep spectral feature, covering a range of energies 
between ~ 0.1-10 PeV x (m^'/MeV)^. Current IceCube 
data already excludes this case for a Z' lighter than a few 
MeV. 

We also emphasize that a Z' can alter the ratios of 
the neutrino flavors that reach Earth, leading to a dif¬ 
ferent distribution of muon tracks, showers, and tau- 
unique events at IceCube. Combining this information 
with measurements of the neutrino spectrum can signif¬ 
icantly extend IceCube’s sensitivity to Z' models and to 
other exotic physics scenarios. 

IceCube’s discovery has opened a new window into the 
interactions of neutrinos at high-energies and over very 
long baselines. As IceCube and other neutrino telescopes 
continue to refine their measurements of this population 
of extragalactic neutrinos, this data will become increas¬ 
ingly sensitive to physics beyond the reach of laboratory 
experiments. As we have shown, IceCube’s current data 
already excludes a small range of the Z' models consid¬ 
ered here. As more data is collected, the range of models 
within the reach of neutrino telescopes will increase, al¬ 
lowing us to explore a significant fraction of remaining 
allowed parameter space. 
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